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Abstract. TauDecay is a library of helicity amplitudes to simulate polarized tau decays, constructed 
in the FeynRules and MadGraph5 framework. Together with the leptonic mode, the decay library 
includes the main hadronic modes, r — > v T + n, 2ir, and 37r, which are introduced as effective vertices by 
using FeynRules. The model file allows us to simulate tau decays when the on-shell tau production is 
kinematically forbidden. We also demonstrate that all possible correlations among the decay products of 
pair-produced taus through a Z boson and a scalar/pseudoscalar Higgs boson are produced automatically. 
The program has been tested carefully by making use of the standard tau decay library Tauola. 



1 Introduction 

The recent observation of a standard-model-like Higgs bo- 
son with mass around 125 GeV at the LHC pQ has in- 
creased our interests in the H — > t + t~ decay mode, which 
has the modest branching fraction yet to be established [5] • 
On the other hand, tau leptons have been well known as 
a particle which can explore the Higgs sector [3113] as well 
as supersymmetric models through their polariza- 

tion [IBS]. 

Tauola [TOIimiT2"] is a well-known program to simu- 
late polarized r decays, including various hadronic modes, 
and recent theoretical and experimental event simulations 
involving r decays heavily depend on itQ The standard 
Tauola Universal Interface [IB] includes the longi- 
tudinal spin correlations between pair-produced taus [T7] , 
while the extended version also takes into account the 
transverse spin effects [T8l[T9] . We point here some lim- 
itations of the use of Tauola: 

i) The produced r has to be stable, i.e. on-shell, in event 
generators. 

ii) The interface for the spin correlation effects is limited 
to the standard processes such as Z/j* — » t + t~ and 
H/A^t+t-. 

a e-mail: tong.li@monash.edu 

b e-mail: kentarou.mawatari@vub.ac.be 

c e-mail: junnaka@post.kek.jp 

1 The event generators, Hekwig++ and Pythia8, have their 
own r-decay package [1311 14] : TauSpinner for studies on spin 
effect in tau production was recently reported [15] . 



hi) For the transverse spin effects, a particular simulation 
is required by accepting or rejecting a pair of produced 
taus based on a spin weighting factor. 

In this article, we present a new implementation of a 
r-decay model file, and construct a library of helicity am- 
plitudes, TauDecayJ^] to simulate polarized r decays in 
the FeynRules (FR) 20 and MadGraph5 (MG5) [21] 
framework. In Section [2] we describe how we implement 
the effective vertices via FR for the main hadronic de- 
cay modes, namely r — > v T + tt, 2ir, and 37r, and present 
the validation of our TauDecay package. To address the 
limitations mentioned above, in Sect. [3] we study scalar- 
tau (stau) decays in scenarios with stau being nearly de- 
generate in mass with neutralino, where the on-shell tau 
production is kinematically forbidden. Moreover, in Sec- 
tions [4] and [5] we demonstrate that all possible correlations 
among the decay products of pair-produced taus via a Z 
boson and a scalar/pseudoscalar Higgs boson can be pro- 
duced within our full-fledged package. Sect. [6] presents our 
brief summary. 

2 Effective vertices for hadronic tau decays 

In this section, we consider effective vertices for the three 
hadronic r decay modes, which together with the leptonic 
mode account for 90% of the r decays [55], and describe 
how we implemented them via the FeynRules (FR) [20] 
and MadGraph5 (MG5) [ST] packages. 

2 The TauDecay package is available at 
http : //madgraph . kek . jp/KEK/ 
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Gf Fermi constant 1.166379 x 10" 5 [GeV -2 ] 
f-n pion decay constant 0.13041 [GeV] 
cos 8 C Cabibbo angle 0.97418 

Table 1. Physical constants 22 . 



The decay modes which we consider in this article arc 



7t mode : r~ 
p mode : r~ 
ai mode : t~ 



is t tt 
v T p~ 



V T TT 7T°, 
lS T lT p~ 



(la) 
(lb) 



^ r 7r 7r°7r", (lc) 
r~ — > v T a\ — > v T it~ p° — > V T TT~ 7T~ 7T + , (Id) 

as the 2ir and 3n modes are dominated by the p and a\ 
vector-meson productions, respectively. 

The t-v t -ti vertex can be introduced by the effective 
interaction Lagrangian: 



£ w = V2G F ,h rfP L v r 9u7r- + h.c 



(2) 



with the physical constants in Table [T] the chiral projec- 
tion operator Pj,, and the constant form factor 



fl = U cos9 C - 
The effective Lagrangian for the p mode is 

C p = V2G F f 2 f^P L V T (TT^TT- - TT-^TT ) + h.C 

with 

f 2 = V2cosd c F p (Q 2 ). 
The form factor F p {Q 2 ) is parametrized by [2"3"IH"U] 

F p (Q 2 ) = [B p (Q 2 ) + aB p ,(Q 2 )]/(l + a) 
with the Breit-Wigner factor 

2 

„ mi. 

B V (Q 2 ) = 



rr v (Q 2 )' 



(3) 
(4) 

(5) 
(6) 

(<) 



where Q = qi+ q 2 for r 
width is 



r v {Q 2 



vit {qi)^ (jq 2 ) . The running 

y/Q* gv(Q 2 ) 



my g v (jn 2 v y 
where the p meson line shape factor is 



9 P 



(Q 2 ) = ?( 



Q 2 ' Q 2 



with 



(3(a, b) = (1 + a 2 + b 2 - 2a - 2b - 2ab) 1/2 . 
We take a = -0.145 in ^ [23]. 



(8) 



(9) 



(10) 





Fig. 1. Effective vertices generated by MG5. 

In practice, we introduce only pions, ^ and 7r°, as 
new particles and implement the above Lagrangians ^ 
and Q into FR^] which provides the Ufo (Universal 
Feyn Rules Output) model file [21], with f 2 as a con- 
stant parameter. The Aloha (Automatic Libraries 
Of Helicity Amplitudes) program [25] in MG5 reads 
the model file to create the Helas (HELicity Ampli- 
tude Subroutines) library [25] for helicity amplitude 
computations. At this stage we replace the constant pa- 
rameter by the momentum-dependent form factor |5]). 

Unlike the above two cases, the effective vertex for the 
a\ mode cannot be obtained by the Lagrangian since the 
vertex structure is not symmetric between the two identi- 
cal pion in the final state. Therefore, instead of introduc- 
ing the Lagrangian, we implement it by hand in FR and 
created the Ufo model file|^] The effective vertex, or the 
decay amplitude, for the a\ mode is 

M ai = V2G F f^P L v T J», (11) 
where the hadronic current J M is given by [2"3lfTU] 

^ = f 3 [F 13 (q^ - qg - G 13 Qn + 2)] (12) 



with Q — qi + q 2 + q 3 for r 
The form factors are 



v t tx (qi)n {q 2 )ir+(q 3 ). 



f 3 = cosOcB^Q 2 ), 



(13) 



F P (Q 



g i3 = Q-fe- g 3) ; (14) 



with Qa = qi + q 3 (i = 1,2). The a\ meson line shape 
factor in (l8|) is parametrized as [231H0] 



2^3 



9a 1 (Q 2 )={ 



4.1 

Q 2 

x [1 - 3.3(Q 2 - 9ml ) + 5.8(g 2 - 9m 2 ) 2 } 

if Q 2 < (m p + m^) 2 , 
10.38 9.32 0.65 

if Q 2 > (m p + m n ) 2 . 

(15) 



1.623 



Q 2 



The 7r°7r°7r _ mode is the same, except the mass difference 
between ir^ and ir°. 

Our Ufo model now allows MG5 to generate Feyn- 
man diagrams including hadronic r decays (Fig. [T]) and 

3 7T° decay is not considered in the program. 

4 Input option of the function WriteUFO in FR allows us to 
include such a vertex [24] . 
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particle 


mass [GeV] 


width [GeV] 


T 


1.77682 


2.265xl0~ 12 




0.1395702 




7T° 


0.1349766 




P 


0.77549 


0.1491 


p(1450) = p 


1.465 


0.4 


ai 


1.23 


0.42 



Table 2. Particle masses and widths 1221 



mode width [10 _i;i GeV] 

PDG [221 Tauola [12] TauDecay 



e~v 


4.04 


4.02 


4.04 


pL~ V 


3.94 


3.91 


3.94 




2.45 


2.47 


2.42 


TV-n 


5.78 


5.39 


5.39 


7r°7r°7r 


2.11 


2.25 


2.27 




2.04 


2.21 


2.22 



Table 3. r decay partial widths. 
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Fig. 2. 7r 7T° (left) and 7r°7r°7r (right) invariant mass distri- 
butions for the t decay into two and three pions. 



the corresponding helicity amplitudes. TauDecay is a li- 
brary to simulate polarized t decays, based on the r decay 
helicity amplitudes created by MG5 with the form factor 
implementation. Using TauDecay, we produce the par- 
tial decay widths in Tabl e [3] as well as the pion invariant 
mass distributions in Fig. |2j The results from the standard 
r decay library Tauola ^12] are also displayed for com- 
parison. We note that the difference in the large m n -„o 
invariant mass region for the p mode comes from the differ- 
ent parameter choice for m p ', which is taken as 1.37 GeV 
in Tauola. Moreover, the fractional energy distributions 
of the polarized r decays in Fig. [3] agree with the ones in 
the collinear limit [3] as well as by TAUOLA. 



3 Stau decays via an off-shell tau 

An important application of our FR r-decay model file 
is for the case that on-shell r production is kinematically 
forbidden. The model file allows MG5 to treat an inter- 
mediate r as a propagator. 

Such a case can happen in the so-called stau-neutralino 
(fi-Xi) coannihilation scenario in the constrained mini- 
mal supersymmetric standard model (CMSSM) . The neu- 
tralino x? is the lightest supersymmetric particle (LSP) 
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Fig. 3. The fractional energy distributions of p and ai (top) 
and 1-prong ir~ (bottom) from left-handed r~ (left) and right- 
handed r~ (right) at E T — 50 GeV, normalized to the respec- 
tive branching ratios. The leptonic decay mode is also shown. 
The energy fractions are z = E p , ai /E T (top) and E n , e /E T 
(bottom) in the laboratory frame. 



and stable by assuming the i?-parity conservation. The 
scalar-tau (stau) fi is the next-to-lightest SUSY particle 
(NLSP), which exclusively decays into a r-lepton and a 
LSP Xi- The scenario is cosmologically preferred to pro- 
vide the observed dark matter relic density [57] as well 
as to solve the 7 Li problem in the standard big-bang nu- 
cleosynthesis [28]. The recent global fit analysis for the 
CMSSM also favours the point [5^ . If mf 1 — m%o > m T , 
there is no obstacle to use Tauola for r decays, following 
event generation with r-leptons as an on-shell particle [SU] ■ 
For the case of 



Am 



A.1 



(16) 



however, we cannot generate events unless r decays are 
taken into account in event generators. One of the common 
benchmark points for the SUSY searches at the LHC is e.g. 
CMSSM40.1.2 [3T], which provides the above mass spec- 
trum with m fl = 230.3 GeV and m^o = 228.7 GeV [32] Q 

When the two-body decay mode f x — > Xi + r is closed, 
the three-body or four-body decays via the off-shell r (see, 
e.g. Fig. Q become dominant and the fi could be long- 
lived |35ll36ll2"9"] . Figure [5] shows the stau partial decay 
widths for each r decay mode against the x? mass. The fi 
mass and the relevant mixing angles in the neutralino and 
stau sectors are fixed as in CMSSM40.1.2, providing a f^- 
like stau and a bino-like neutralino. The vertical dashed 
lines denote the threshold for each decay mode, where 
a\ and p indicate Am = m^ and m2n, respectively. As 
the Xi mass is approaching the t\ mass, i.e. Am becomes 
smaller, the decay width becomes smaller very rapidly. 
Just below the r threshold, around the CMSSM40.1.2 



5 A similar situation in the chargino-neutralino degenerate 
scenario has been studied in 33 : 34 . 
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Fig. 4. Feynman diagram of a four-body fi decay generated 
by MG5. 
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Fig. 5. Stau partial decay widths as a function of the LSP 
neutralino mass. The sum of the partial widths is also shown by 
a black dashed line. The stau mass is fixed at 230.3 GeV, based 
on the CMSSM40.1.2 benchmark point. The vertical dashed 
lines denote the threshold for each decay mode. 



point (denoted by a vertical solid line), decay widths of 
all the modes are comparable with those in the on-shell 
r case, while the 7r mode is dominant in the small Am 
region due to the phase space suppression for the multi- 
pion and leptonic modes. Only in the Am < m^ region 
the electronic mode becomes dominant. Those agree with 
the recent study [29j , although p and a\ are considered as 
on-shell particles in their calculations. 

The collider signature of the long-lived ft significantly 
depends on the f\ width, i.e. the lifetime; sec e.g. Refs. 37, 
[551155] , and as shown in Fig. [5] the lifetime strongly de- 
pends on Am. For instance, the f\ predicted at CMSSM40.1 
whose lifetime is O(10~ 8 s), could leave a charged track 
with a displaced vertex of its decay inside the detectors. 
We briefly study the decay distributions for such a long- 
lived stau to see the left-right mixing of ft in this situation. 

The scalar partners of left-handed and right-handed 
taus, fi and tr, mix to form two mass eigenstates, and 
the lighter one is 










p 








J 







0.0 0.5 1.0 1.5 2.0 

E[ GeV| 

Fig. 6. Energy distributions of the one-prong pion in the 7r 
mode (left) and in the p mode (right) for the left-handed f 
(blue) and right-handed f (red) decay in the stau rest frame. 
The distributions for the intermediate p is also shown by 
dashed lines as a reference. 



where 6f is the mixing angle. Similarly gauginos, B and 
W3, and neutral Higgsinos, H% and H®, mix to form four 
mass eigenstates, and the lightest one is 



(18) 



where X = (B,W 3 ,H%,H°). The chirality of r in the 
f\ — > TXi decay depends on these mixings and determined 
by the interaction Lagrangian 



C = xi{a L P L + a R P R )Tf* + h.c. 



(19) 



n 



cos Uf tl 



TR, 



(17) 



with the chiral-projection operators, Pr/l — ^(1±7 5 ). 
The couplings, cll and cir, are 

= cos9f-^=(U2i + Un tan8w) + sm8f UziY T , (20a) 
v2 

clr = sin6?~2U n tan8 w - cos6» f U 31 Y T , (20b) 
v 2 

where Y T — — gm T j\[2m w cos/3, g is the SU(2) L gauge 
coupling, and tan j3 is the ratio of the vacuum expectation 
values of the two Higgs doublets. For most of SUSY sce- 
narios as well as for CMSSM40.1.2 the lightest neutralino 
is gaugino-like (C/31 ~ 0), where the chirality is conserved. 

Figure 6 shows the energy distributions of the one- 
prong pion in the it mode (left) and in the p mode (right) 
for the left-handed f (blue) and right-handed f (red) de- 
cay in the stau rest frame. The distributions for the re- 
constructed p is also shown by dashed lines as a reference. 
The phase space of the decaying off-shell t is limited and it 
is no longer energetic. Therefore the energy distributions 
2of the pion are quite different from the ones in Fig. [3j 
where the collinear limit is safely applied, although the 
remnant of the difference between left-handed and right- 
handed can be still seen especially for the n and p modes. 



4 Spin correlations in tau pair decays 

Another important application of the FR r-decay model 
is spin correlations in tau-pair production BMH 
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IT3 | fT4| , fT5] jj As in the previous section, one can generate 
amplitudes including the r decays, and hence the spin 
correlations between the two taus are automatically taken 
into account. 

In this section, we present correlations in tau-pair pro- 
duction through a Z boson, a CP-even scalar H, or a 
CP-odd scalar A, 



qq 
99 
99 



H 



(21a) 
(21b) 
(21c) 



and the subsequent hadronic r decays. We note that, al- 
though we consider the above standard r-pair production 
processes in this paper, our r-decay program is indepen- 
dent of the t production and can cooperate with all kinds 
of new physics models. 



4.1 Helicity formalism 

First, we discuss the case that both taus decay to v T ir by 
using analytic expressions. We assign the four-momentum 
and the helicity of each particle for the process as 

ai(pi,o-i) + a 2 {p2,<y 2 ) 

t~((7i,Ai) + T + (g 2 ,A 2 ) 
->■ Tr-(fci) +ir+(k 2 ) +v r (ka,-) + v T (k 4 ,+), (22) 

where ai.2 stand for quarks or gluons (or even electrons 
or photons in case of e + e~ or 77 collisions). The helicities 
take the values 0^/2 and A,/2 for quarks and leptons and 
Oi for gluons and photons with = ±1 and A.; = ±1. 
Although we evaluate the full amplitude for the above 
processes numerically to present the distributions, we dis- 
cuss it below in terms of the r-pair production amplitude 
and the r decay amplitudes, which give us better under- 
standing of the distributions. Moreover, we construct our 
TauDecay package based on the following analytic ex- 
pressions in the end of this section. 
Using the completeness relations 

ji + m = y]u(qi, Xi)u(qx, Ai), (23a) 

Ai 

i2-m = 2_ J v[q 2 ,\2)v(q2,\2), (23b) 

the full amplitude can be expressed as the product of the 
tau-pair production amplitude (Mx=z,h,a) and two r — > 
ttv decay amplitudes (.M1.2): 

M(pi,ai;p2,a- 2 ;ki) 

= D {qi) D ((i) E ■ M x(pi,<J 1 ;p 2 ,(T2;qi, Ai; g 2 , A 2 ) 
x Mi(qi,Xi;ki;k 3 ) M 2 (q2, A 2 ;fc 2 ;fc 4 ) 



(24) 



6 The correlation of r r from a doubly charged Higgs bo- 
son H has also been studied in 1401. 



with the r propagator factor D(q 2 



i 2 + im,r) 1 . 



It is straightforward to obtain the squared matrix ele- 
ments of the full production plus decay amplitudes, 



Y,\M\ 2 = ^2\M(p 1 ,a 1 ;p 2 , ( r 2 ;k i )\ 2 

Cl,2 

= i^m<72 2 )| 2 e E^t^; 2 ^ 



Ai,A2 Ai,A2 



(25) 



in terms of the production density matrix V^ 1 ^ 2 and the 



decay density matrices Vi^ 



Ai 2 



Ai A2 



E-M^ 2 (Ml 1 ' 



V^=M Xl (M- Xl Y 
V 2 ^=M X2 (M- X2 Y 



(26) 

(27) 
(28) 



In the narrow width limit the propagator factor becomes 



\D(q 2 )\ 2 ^-^S(q 2 -m 2 ). 



(29) 



Although we only consider parton-level subprocesses, one 
can generalize (25 1 to mixed case and apply it for any pro- 



cesses, including summation over different subprocesses 
and a product of the relevant parton densities. We can 
also easily replace the r — > v T i: decay amplitude by one 
for the other hadronic decay mode as well as the leptonic 
mode. 



4.2 Kinematics 

Let us define the kinematical variables. In the collision 
center-of-mass (CM) frame, i.e. in the X rest frame, we 
choose the r momentum direction as the z-axis, 

pi = ^(1,- sin 0,0, cos 0), 
p 2 = ^(1, sin 0,0,- cos 0), 
^■,0,0,0), 



1\ 



<12 



^(1 + ^,0,0,-/3), 



(30) 



where /3 = f-) with /3(a,6) defined in ([Toj), is 

the scattering angle, and we choose p\ x q\ direction as 
the y-axis. The momenta of the t~ decay products are 
parametrized in the t~ rest frame, 



h = 



: (1 



f3i sin 9\ cos 4>i , /3i sin 6 1 sin </>i , f3i cos 9\ ) 

(31) 



2 V- ■ qt 

fc 3 = ^^-^(l, — sin 6\ cos — sin 0i sin^i, — cos#i), 
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The amplitudes are non-zero only when a = o\ = a 2 
and A = Ai = A2. Here, y T = v2m T /v is the r Yukawa 
coupling, and 175 is the ggS coupling, which is given by 
9h = a s gHtt/3nv and qa = a s gAtt/^v in the heavy-top 
limit. 

The t — y i^v decay amplitudes in the r rest frame are 



Fig. 7. Schematic view of the coordinate system. 

with (3i — 1 — m 2 ,/q 2 . Similarly, those of the r + decay 
products are 



where 



M Xl = a/1 + Ai 1 



cosfc/i e lAl0l/2 , 
M X2 = v/l + A 2 cos^e^ 2 ^ 2 . 



(38) 

(39a) 
(39b) 



^■(1 + -^f,^2 sin6> 2 cos 02, P2 sin 9 2 sin 4> 2 , f3 2 cos 9 2 ), 



"r. 



$2 (1 , — sin 02 cos 02 , — sin 6> 2 sin <fi 2 , — cos 9 2 ) 



(32) 



with /3 2 = 1 — rn 2 /q 2 . We note that the polar (2-) axis 
and the y-axis normal to the scattering plane are chosen 
common to all the three frames, and the two decay frames 
differ only by the boost along the r production axis (see 
Fig. Fr}. The r width is very narrow, r ~ 0(lCr 12 GeV), 
and hence we take the narrow width limit, q\ — q 2 = to 2 , 
in the following analytic amplitudes. 



4.3 Tau-pair production and decay amplitudes 

The helicity amplitude for the tau-pair production via Z- 
boson (21a) is given by 



M 



e 2 s 



, : - M XlX2 



s ~ to| + imzTz 
where the reduced amplitude is 
M XlX2 

-\g%{g T +{l + XiP) + g T _(l - \ip)}(<r\i + cos0) 

for Ai 7^ A 2 , 
for Ai = X 2 - 



(33) 



^X ig i(gl+g-) S m0 



(34) 



Here, we neglect the initial fermion mass and take a — 
o\ = —o 2 . g q ± and g T ± are the Z-boson couplings to right - 
and left-handed quarks and tau leptons, respect ively . 
The amplitude for the S(= H, ^4)-boson case ( 21b ) and 
is 



(21c 



M s = 



-Vrgss 



s — m s + imsTg 2 



M 



\i_\-2 



(35) 



where 



M xx = (3X 
M xx = i 



for S = H, (36) 
for S = A. (37) 



Those amplitudes are invariant with the boost to the col- 
lision CM frame. 



4.4 Helicity correlations 

The trivial helicity correlations are given by the diagonal 
parts of the density matrices, Ai = Ai and A2 = A2. In 
the following analytic expressions we take m T /\ps = 0, 
i.e. [3=1. 

For the Z production, only for Ai 7^ A2 the production 
amplitude is non-zero. After integrating out the scattering 
angle and the azimuthal angles 01,2, the squared matrix 
clement (25) is 



n --n , VZX^>\-V 2 ~ 

oc (1 — cos 61 cos 9 2 ) + k (cos 61 — cos 9 2 ). (40) 



vxzv x xv 2 



withK = {g T + 2 -gl 2 )/{g T + 2 +g'/ 2 ) 0.15. As seen in 

the polar angle cos#i j2 distribution of the pion arising 
from is (1 — cos #1,2), while from r^jr it is (1 + cos#i j2 ). 
Hence in Z — > t^t^ or t^t^ decays the two pions tend 
to be emitted to the opposite direction along the z-axis. 
The difference between the left and right coupling, i.e. 
parity violation, gives a small linear dependence of cos 6\ 
and cos 9 2 in the single differential cross section, and one 
can see a slightly higher density in the cos6*i < and 
cos 9 2 > region in the left panel of Fig. [8] 



H/A 



■ 

'.V. :■.•>■■■:■'., ■ ..-'r/.. 



1 « m 



1 ■ 



cose, 



Fig. 8. cos(9i-cos#2 correlation in pp — > X — > r (— > 
ir~v) r+(^ tt+P) for X = Z (left) and H/A (right). The pion 
polar angles cos 61,2 are defined in the r rest frame; see Fig. [7] 
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For the Higgs production, on the other hand, only for 
Xi = X 2 the production amplitude is non-zero, and there- 
fore 



VZZT>xZT> 2 Z oc (1 + cos0icos0 2 )- (41) 



Apart from the linear dependence of cos #12 for the Z case, 
completely the opposite is favoured for H/A — > t/tZ or 
t r t r d eca y s as shown in the right panel of Fig. ^1 Note 
that the helicity correlations for H and A are identical. 
For comparison the masses of H and A are assumed to be 
the Z-boson mass, m^/A = 



4.5 Polarization correlations 

The non-trivial polarization correlations are given by the 
off-diagonal_parts of the density matrices, i.e. Ai = — Ai 
and A2 = — A2, which produce the azimuthal angle depen- 
dence. When we isolate the azimuthal angle dependence in 
( 25 ) , there are nine distributions (including one constant 
piece) as 



pAiA 2 v x 1T> \ 2 

A1A2 Ai X-2 



F 1 + {2Ke[F 2 cos 
+ (5Re 



t + F 3 COS 02 

F± cos(0i + (j> 2 ) + FZ cos(0i - 



<h)] 



sm 



)}• 



(42) 



Here, and in the following, summation over repeated in- 
dices (Ai, A 2 , Ai, X 2 ) = ± is implied. The coefficients 
are the functions of the kinematical variables except the 
azimuthal angles <fii y2 . F° r the production of a Z boson, 
they also depend on the production angle 0. For the spin- 
particle case, only the FZ term in ( 42 ) survives due to 
the helicity selection Ai = A2. All the sine terms vanish 
when CP is conserved and when the absorptive part of 
the amplitudes are neglected, e.g., in the tree-level ap- 
proximation. 

Because the phase of the product of the two decay 
density matrices is 

2?i£x> 2 £ oc exp[i{(A! - Ai)& - (A 2 - A 2 )<M/2], (43) 

the coefficients F^_\ are expressed in terms of the produc- 
tion density matrix and two decay density matrices as 

A 2 



F 2 =VZ%V 1 ±V 2 %, 
F ± 



(44) 



The azimuthal angle correlations are manifestly expressed 
by quantum interference among different helicity states of 
the intermediate tau leptons. 

The Z production could produce the <j>\ (or (f> 2 ) depen- 
dence, however this is very small because F 2 (P3) is pro- 
portional to m/y/s. Similarly, F± oc 0(m 2 /s), and hence 



r'dr/d* 






H A „.,,.»» 




z 




Fig. 9. Azimuthal angle correlations, <j)i + (f>2 (left) and cf>i — 4>2 
(right), in pp — > X — > t~(— > h~v)t + (— > -k + v) for X — Z 
(blue-dotted), H (black-solid), and A (red-dashed). The pion 
azimuthal angles 0i,2 are defined in the r rest frame by the 
pp — > t + t~ scattering plane; see Fig. [7] 



the distribution of the azimuthal difference between the 
two r-decay planes, 4>i — <f> 2 (= A<fr), is flat. An interesting 
correlation between the two decay planes for the Z case is 
the 4>i + 4> 2 (= #) correlation, whose coefficient F± is 



Fa 



-g\g T _ sin 2 sin 9\ sin 9 2 . 



(45) 



After integrating out 0, 6\ and 9 2 , the azimuthal asym- 
metry is given by 



with 



1 dr 



9 F + 
4 F 1 



1 

2tt 



[l + Aj cos<2>] 



-9+9- 
2 +9-'' 



0.30. 



(46) 



(47) 



The distribution is enhanced around <P — and 27r, while 
it is suppressed around <P = tt. 

For the H/A production, as mentioned above, only the 
<pi — (f> 2 term is non-zero and the coefficient is 



F 4 cx qp sin 0\ sin 6 2 , 



(48) 



where the — /+ sign is for CP-even/odd scalar. The az- 
imuthal asymmetry is given by [3] 



1 dr 1 r , „ n 



with 



A] 



2FZ_ 
Fx 



16 



T0.62. 



(49) 



(50) 



The CP-even and -odd scalars have opposite modulation, 
and the distribution is strongly enhanced (suppressed) 
around A<j> = tt for the CP-even (-odd) scalars. 

To examine the validity of the model file, Fig. [9] shows 
the normalized azimuthal correlations between tt~ and 
7r + , and our numerical results agree well with the above 
analytic formula. 
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Fig. 10. The fractional energy correlations in the laboratory frame in pp — > X — > r (— > 7r i/) r + (^ A + v) at 7s = 8 TeV for 
X = Z (top) and H/A (bottom), where A + is the leptonic, ir, p, and ai decay modes from the left to right. The resonant mass 
is taken to be the Z-boson mass. 



4.6 Spin correlations at the LHC 

To present the validation of our program, we have so far 
discussed the angular distributions, assuming that all the 
T-decay products can be reconstructed. Here, we show 
some distributions of the realistic experimental observ- 
ables at the LHC. 

First, we extend the helicity correlation in r — » i/ T 7r 
discussed above to other tau decay modes. In Fig. [10] 
we present the helicity correlations in terms of energy 
fractions in the laboratory frame for pp — > Z/H/A — > 
T~(—> -k~v)t + (-^ A + D), where A + is the leptonic, tt, 
p, and a,\ decay modes, respectively, from left to right. 
The one-prong r decays are considered and the energy 
fractions are z\ = E n -/E T - and z-i = E^ + /E T + with 
A = e, 7T, n(p), 7r(ai). For comparison, the resonant mass 
is taken to the Z-boson mass, and hence the both taus are 
highly boosted, where the collinear approximation can be 
safely applied. For the A — tt mode, in the collinear limit 
(f3 = 1) we have 



r _1 dr/ dcoss 




Higgs rest frame 




A / r 


h L L 







r 'dT/dcosS 




lab frame 






FrfU-Jl H 




a ! -L: 



0.94 0.95 0.96 0.97 0.98 0.99 1 

5/ji 



0.94 0.95 0.96 0.97 0.98 0.99 1 



Fig. 11. The -k + tt acollinearity distribution in the Higgs 
rest frame (left) and in the laboratory frame (right) for pp — > 
H/A -> r~(-¥ tt~v)t + (^ tt + p) with m H/A = 125 GeV at 
7i = 8 TeV. 



zi = (l + cos0i)/2 and z 2 = (1 - cos0 2 )/2. (51) Fig 



As a result, the z\-zi correlation behaves opposite to cos Q\- 
cos 82 correlation as shown in Fig. [8] The helicity correla- 
tions are different between the spin-1 and spin-0 bosons, 
while those are identical between the CP-even and -odd 
scalars. 

Second, an useful observable to measure the parity of 
Higgs boson is the acollinearity angle, namely the angle 8 
between tt + and tt~ [4|. The azimuthal asymmetry Acj) 



can be given as a function of angle 6. To display the 
measurable difference of Higgs parities at the LHC, in 
we show the 7r + 7r" 
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acollinearity distribution in 
the Higgs rest frame (left) and in the lab frame (right) for 
pp -> H/A -> r _ (-> tt-v)t+(-^ tt+D) at s/s = 8 TeV. 
Here we take the Higgs mass as 125 GeV. One can see the 
distinguishable difference between CP-even and CP-odd 
scalars around 5 = tt, although in the lab frame the boost 
effect leads to relatively unclear discrepancy. While the re- 
construction of the Higgs rest frame at hadron colliers is 
difficult, it is possible in the Z-associated Higgs produc- 
tion at e + e~ colliders and the acollinearity distribution 
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in the reconstructed H/A rest frame has been discussed 

in rrsi. 



5 Spin correlations with TauDecay 

To validate our FR r-decay model file, all the results 
shown in Sections [3] and [4] are generated in MG5 with in- 
termediate taus being as propagators, denoted by method A. 
Unless intermediate taus are off-shell such as the case in 
Sect. [3j a library which collects all the possible r decay 
channels would be much more efficient for practical event 
generations since t production and its decay can be sim- 
ulated independently as we use Tauola. 

TauDecay is a Fortran library of r-decay helicity 
amplitudes, constructed in the FR/MG5 framework as de- 
scribed in Sect. [2] The input of the subroutines for each 
decay channel is the four momenta of the decaying r as 
well as its decay products and their helicities. The out- 
put is the amplitude, i.e. a complex number, at a given 
phase space point with a given helicity configuration. In 
this section, we present how the TauDecay package can 
reproduce spin correlation discussed in Sect. [4} 

We show two methods by means of TauDecay; with- 
out and with spin density matrix, denoted by method B 
and C, respectively. The both methods are carefully checked 
by method A. 

B: [without density matrix) The r production amplitude 
as well as the decay amplitudes from TauDecay are 
evaluated for a given process for a given phase space 
point. The r helicities in the product of the production 
and deca y am plitudes should be summed over as indi- 
cated in ( |24[ ), and then the r helicity summed ampli- 
tudes should be squared. This is theoretically identical 
to the propagator method in the narrow width limit, 
but it is useful for systematically studying different tau 
decay modes for each production process. 




Fig. 12. Azimuthal angle A<fi correlation in pp — > H/A — > 
t~(— >• ix~v) t + {— > -k + v) by three different methods; A. propa- 
gator (solid), B. without density matrix (dashed), and C. with 
density matrix (dotted). 



C: (with density matrix) First, the r production density 
matrix is computed for each phase space point. If there 
is only one tau in the final state, this will be a 2 x 2 
matrix, while it will be a 4 x 4 double density matrix 



for two taus as •p^ 2 in (|26|). Then, the correlated tau 

Ai A 2 1 1 



decays are simulated by multiplying the production 
density matrix with the corresponding r decay helicity 
amplitudes and their complex conjugates provided by 
TauDecay, namely the decay density matrices 2?i^% 



according to (25 1 



Figure [12] shows the azimuthal angle A<j> correlation in 
pp —> H/A —> t~(—> 7r~i/) t + (-^ n + D) for the comparison 
of the three methods, and the same results are reproduced 
as expected. 



6 Summary 

In this paper, we have implemented the main hadronic 
tau decays, t — » v T +n, 2ir, and 37r, by using the effective 
vertices in the FeynRules and MadGraph5 framework, 
and constructed TauDecay, a library of helicity ampli- 
tudes to simulate polarized tau decays. The model file 
allows us to simulate tau decays when the on-shell tau 
production is kinematically forbidden, and the stau decay 
in the stau-neutralino coannihilation region was discussed 
as an application. We also demonstrated that all possible 
correlations among the decay products of pair-produced 
taus through a Z-boson and a scalar/pseudoscalar Higgs 
boson can be produced within our full-fledged package. 
The program has been tested carefully by making use of 
the standard tau decay library Tauola. 
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